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PACS 67. 10. Jn - Quantum fluids: Transport properties and hydrodynamics 
PACS 05 . 40 . Jc - Brownian motion 

PACS 71 . 38 . -k - Polarons and electron-phonon interactions 

Abstract - Recent experiments with cold atoms on the impurity motion in one-dimensional liquids 
of interacting bosons have revealed an interesting interplay between the polaronic impurity mass 
shift and the renormalization of the optical potential. We show that the influence of the external 
trap on the Bose gas leads to a steeper effective potential for the impurity. We propose a framework 
in which this potential renormalization and the mass shift can be quantitatively understood by 
combining a semi-classical theory of density wave excitations in the Luttinger liquid with the non 
equilibrium formalism of a quantum Brownian particle. The obtained theoretical results reproduce 
well recent experimental data. 



Diffusion in low dimensional quantum liquids has been 
a major field of research in the last decade [TUlOj. In onc- 
dimensional (ID) liquids of interacting bosons a moving 
impurity is subject to a drag force [6] and dissipates en- 
ergy for all velocities even at zero temperature [5]. The 
experimental design of artificial ID impurity-quantum liq- 
uid systems has now become possible by confining cold 
atoms in optical nanotubcs. Using these techniques, the 
diffusion of impurity atoms in contact with a Luttinger 
liquid (LL) with tunable impurity-LL interaction was re- 
cently studied |10| . Due to the external trapping potential 
the minority atoms undergo damped oscillations which di- 
rectly confirm that dissipation takes place in this system. 

In Ref. [TT] the authors attempted at a theoretical de- 
scription of this experiment with the Gross-Pitaevskii for- 
malism. In this Letter we pursue an alternative approach 
by considering the impurity atom as a quantum Brown- 
ian particle. The quantum liquid then plays the role of 
an exotic quantum bath and it can be dealt with using 
Luttinger theory. As already pointed out in [TDHT^ the 
impurity atom acquires an effective mass due to its in- 
teraction with the LL. Moreover, the external trapping 
potential leads to an inhomogcncous density profile of the 
LL with non-trivial effects on the impurity motion. In this 
Letter we combine several independent ideas that will al- 
low us to: (a) estimate the mass shift of the impurity, (b) 
evaluate the effect of the non-homogeneous density profile 



of the LL as. to a first approximation, the renormalization 
of the confining potential, (c) use the non-equilibrium for- 
malism of quantum Brownian motion developed in |12j to 
reproduce the data in |10j quantitatively. 

The free Hamiltonian of an impurity with mass M/, 
in an optical trap modeled by a harmonic potential with 
spring constant k, reads 



■Hi 



2Mi 



(1) 



with p and q the momentum and position operators. 

The impurity interacts with an LL which is confined by 
the same potential [10]. For the sake of simplicity we will 
incorporate the effects of the trap on the LL later. The low 
energy excitations of an unconfined ID quantum liquid are 
described by the Tomonaga-Luttinger Hamiltonian 



dx 
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with the two canonically conjugate bosonic fields n(x) and 
(^(x) [9lfT3] . 4> is related to the LL particle density through 
p{x) = po — (l/7r)0'(.T). In the second quantization lan- 
guage the Hamiltonian can be equally expressed in terms 
of bosonic operators b\, and b^ which we define below. The 
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dimensionless coefficient K and tlie sound velocity u to- 
tally characterize the low energy properties of such a ID 
system. For translationally invariant cases they only de- 
pend on the Lieb-Linigcr parameter 7 = Mj^WL/hpQ, with 
Mi the mass of the bosons, fiwi, the strength of the inter- 
action and pq the density of the LL [T^. 

We model the impurity-LL interaction through 



n 



IL 



dxdy U{x — y)p{y)S{x — q) 



(3) 



with U (x) the interaction potential, p{x) the LL density 
and q the impurity position operator. In Fourier space 
we define (f>{x) = L~^^^ J2k s~^''^4>k with k = 2Trn/L and 
n e Z. The full Hamiltonian is then H = Hl + 'Hil +^1, 
with Hi defined in eq. ([1]) and 
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Note that we rescalcd the fields according to (f>k 

{■KK/h)tpk and fffe i-> ^/(KpKK)tlk- In terms of b\, bk 
the (rescaled) field reads 0^ ~ ^/i/2|fc|(&Lfc + bk)- We 
choose Uk = hw/VZ g— "|fc|/2'-'c -ftrith some cutoff wave 
vector ujc/u depending on the microscopic properties of 
the coupling. Equation ([5]) considers only the so-called 
forward impurity-LL scattering. The backward scatter- 
ing potential C/back is omitted since it renormalizes to zero 
for an isolated fixed immobile impurity according to the 
equation 9rC^back = (1 ^ -f^)f^back El- In the following we 
assume that this is also true for a mobile impurity and we 
focus on the case where K > 1. 

In |10| the impurities are initially localized at the center 
of the potential tubes with a laser blade that creates an- 
other harmonic potential well with spring constant kq > k. 
After their subsequent release they undergo stochastic dy- 
namics that resemble the ones of a damped harmonic os- 
cillator. Catani et al. measured the equal-time correlation 
function C(t,t) and they drew the following conclusions; 

(I) The oscillation frequency fi/ is virtually not affected 
by the value of the impurity-LL interaction hw. 

(II) Equations ([I])-® resemble the well-known Frohlich 
polaron Hamiltonian that should result in the impurity 
mass renormalization, Af/ 1— > Mj, as a function of the 
interaction hw. Point (I) then indicates that in parallel 
to the mass renormalization the potential spring constant 
should be renormalized as well, k k*, in such a way 
that = \J K* jM^ remained equal to f2/. 

(III) The initial kinetic energy of the impurity can be 
estimated from the high temperature cquipartition theo- 
rem to be ~ 1//3 (note that ?i/3a/ko/M/ ~ 0.1 in [10]) by 
assuming that the impurity has equilibrated with the LL 
before its release. The amplitude after one oscillation 
should therefore scale as ^ 1/-\/k* when neglecting dissi- 
pation such that K*?^ ''^ 1/ P due to energy conservation. 



Furthermore, \/k/k* ~ ^/mJJmJ due to point (I). The 
increase of k* ~ Mj is clearly observed when hw is ramped 
up (see Fig. [1]). Note that for w/wl > 5 the ID regime is 
not ensured any longer which explains the "saturation" of 
K* for large values of w (not described by the effective ID 
theory) . 

(IV) The final (equilibrium) width of the impurity cloud 
is independent of hw. 

We notice that (IV) is at odds with (II). From the 
theory of quantum Brownian motion we know that 
limt_j.oo C(f, t) ~ 1/(/3k*) for a harmonic potential with 
spring constant k* [13]; therefore, the dependence of k* 
on hu should entail a dependence of the cloud width with 
the same parameter. Accordingly, a more thorough analy- 
sis of the coupled systems is needed to correctly interpret 
the experimental evidence. In the following we examine 
the points (I)- (IV) in detail and we propose a way out this 
conundrum. 

The polaronic mass shift. It is well-known that a 
charged particle acquires an effective mass when it inter- 
acts with lattice vibrations through a Coulomb potential. 
Equations dl])-® describe such a polaron with the only 
difference that the interaction is not Coulomb-like. In the 
following we estimate the polaronic mass shift in our prob- 
lem by using the equations of motion (EOM) for cf>k{t) and 
q{t) (see [15] for the use of EOM in this context): 



Mt) + u'k^Mt) = iukd—U^e 
V nh 

Mi^{t) + Kq{t) = 



ikq{t) 



(6) 
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Suppose that the impurity is not accelerated during a 
small time interval, then we can make the Ansatz 

q{t) = m + vt (8) 

0. The solution to 



with V = p{0)/Mi and [q{t) , (j)k(t) 
eq. ([6]) for small t reads 

^k{t) = Ak{t- «)e'^'=** + gfce™'^* + /ifce"™^* 
with the coefficients 

^ikq(0)+ihk^t/2Mj 

Ak{t;v) = ^uk U*J - „ , 



(9) 
(10) 
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0fc(O) + ^0fe(O) - ifc(0; v) - ifc(0; v)- 
ink u 

4>k{Q) - ^4(0) - ifc(0; v) + ifc(0; v)- 
ink u 



The first term on the right-hand-side of ([9]) represents a 
density cloud that moves together with the impurity, while 
the two last terms describe the very wave excitation. For 
instance, in the limit — ^ 00 we have for a mobile impu- 
rity with constant velocity 



/5(x,i) ~^^fc^(^;«)e''^* 



t — ika 



six-m), (11) 
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Fig. 1: ^JkJkF versus the impurity-LL coupling w/wl- Lines 
obtained from \12\ with (from right to left) 7 = 0.25, 7 = 0.35 
and 7 = 0.5. Experimental data points are taken from [10| . 
Triangular points: Result from Feynman's variational theory 

m. 
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Fig. 2: ^C{t,t) for ili/ojc = 2.5, 7 = 0.45 and w/wl = 1(4) 
in the upper (lower) image. Points: Experimental data from 
[To] . Lines: Solution to ((H]) with il} = 



meaning that the LL density profile follows the impurity, 
thus creating a dressed local impurity. It is important to 
note that this simple picture has to be altered when the 
impurity is accelerated. We will come back to this point 
later. 

The combined system of the impurity and the density 
cloud has an effective mass which exceeds the bare im- 
purity mass. As an illustration, we consider the initial 

conditions 4'k{0) ~ Ak{0;Q) and (j>k{0) ~ 0. If the im- 
purity is immobile [i.e. q(t) = q{0)] these initial condi- 
tions lead to the static solution ipkit) = ^A;(0;0) which 
describes a static density cloud without wave excitations. 
Suppose now that the impurity is instantaneously acceler- 
ated to some constant velocity v, then we obtain from ^ 
4>k = Ak{t,v)e'^'^'"* — iAk{0;0)^sinukt. Hence, upon accel- 
eration energy is carried away by a wave excitation, such 
that the kinetic energy of the impurity is less than the ex- 
ternal energy provided. To be more specific, by using ^ 
the average energy of such a wave excitation is found to 
be Ek = ^iv/u)^Ak{0;0)A_k{0]0) with = {v^). We 
define the effective impurity mass through 



M 



J - (l + /i)M/ 



(12) 



with the interaction-dependent correction /i = 
2hw^KuJc/{TT^Mju*). Then, by providing an amount 
of energy E, the impurity acquires after acceleration a 
(mean) velocity given hy E = Mjv'^/2. It is straightfor- 
ward to generalize this calculation to the case where the 
impurity has already a velocity 'Do before the acceleration: 
One simply replaces /i by 



where we used the classical mean value Wq instead of Wq . In 
the following we consider Mj as the true impurity mass. 

The potential renormalization. Ref. [TD] indicates that 
the spring constant of the optical trap is renormalized as 
well. We will show here that the effect of the external 



potential on the LL indeed leads to a renormalization of 
the potential felt by the impurity. In the same spirit as in 
the previous paragraph we study the effects of the external 
potential, which we previously neglected, via its action 
on the density cloud. The force exerted by the external 
harmonic potential on the density cloud is given by 



F : 



uK r- 

= ^ft(^2_^2) >/^^0 , (14) 



where we used ([9]) to express p(t). By considering the 
combined impurity and the density cloud system as one 
entity, F acts in the end on the impurity itself. Interest- 
ingly enough, F changes sign when v exceeds the sound 
speed u such that the subsonic and supersonic regimes are 
quantitatively different. In |10| the impurity moves with 
supersonic speed (\/ (w^) ~ 8.5mm/s while u « 3mm/s 
typically) so that F leads to an increase of the effective 
external potential. In the following we approximate v'^ by 
its mean value = (w^). 

Intuitively, the potential renormalization can be easily 
understood. When the impurity creates a density exci- 
ton it has to push the LL atoms up the optical potential 
to be able to create the density cloud. Therefore it loses 
more energy than what the density wave would cost. The 
inverse is true as well. By absorbing an exciton the im- 
purity gains more energy than the exciton provides since 
potential energy is freed during the absorption process. 
Equation ((T4)) leads to the effective spring constant 



K* ^ {l+fl{v))K 



(15) 



(13) where //(w) = ^^u/ {v^ ^v?'). In conjunction with cq. (1131) 



we thus obtain for the effective potential frequency 



1 



(16) 



In Fig. [T] we compare the prediction in (fTS)) to the ex- 
perimental data [TU]. The best curves arc obtained for 



P-3 



J. Bonartl L. F. Cugliandolol 



7 « 0.25—0.35. Note that all the constants are determined 
by the experimental setup. However, since it is difficult to 
define 7 for a non-homogeneous density we plotted results 
for 7 = 0.5 and 7 = 0.35 and 7 — 0.25 for illustration. If 
the non-homogeneous density profile is approximated by 
an homogeneous one the parameters used in the experi- 
mental setup in [TU] lead to 7 « 0.45. 

Long time dynamics: In the previous two paragraphs 
we studied the potential rcnormalization and the impu- 
rity mass shift by assuming that the impurity velocity was 
constant. Only then is the impurity cloud perfectly local- 
ized around the impurity position. It is clear that such an 
approximation can only hold for short times in the system 
we consider. For instance, while (fTS)) can still be consid- 
ered as a realistic approximation up to the first oscillation 
maximum of the impurity, it certainly fails to describe the 
correct physics for t 00. The general solution to © is 



ds 



smuk(t — s) 
uk 



,ikq{s) 



wave excitations 



(17) 

which leads to a density cloud of the form p{x,t) ~ 
Jp ds S[x — q{s) — u{t — s)] — S[x — q{s) + u{t — s)]. In 
the case of an exponentially damped oscillating impurity 
this density cloud depends only on past values of q when 
t is large and hence, for t — 00, the influence of q{t) on 
p{x,t) becomes negligible. 

To put it in other words, for i — > 00 the density cloud is 
independent of the impurity such that its dynamics decou- 
ple from those of q{t): The LL has no dynamical effects 
on the impurity and one concludes that the LL neither 
renormalizes the impurity mass nor the external poten- 
tial. Accordingly, the final width of the impurity reads 



for t 



(18) 



and not 1//3k*. We have thus found that dynamical quan- 
tities depend on the rcnormalized values k* and Mj while 
final equilibrium quantities have to be computed with the 
bare values k and Mj. We insist on the fact that this 
behaviour has been observed by [lOj where the final impu- 
rity position width is not rcnormalized in contrast to the 
potential rcnormalization that is observed at short times 
(see Fig.[T]). 

The impurity influence functional: We now use the 
Keldysh formalism to derive an effective out of equilibrium 
action for the dynamical impurity position. The action of 
the free oscillator (described by Hi with the parameters 
Mj and k*) is complemented by 
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(20) 



where (Z^(s), <l~{s) are the dynamical Keldysh branches 
with (7^(0) = qi and q~{0) = q'l, and qoir) is the path over 
the initial condition (with imaginary time r) |12j . The last 
line in the right-hand-side of (pil)) describes the initial lo- 
calization due to the laser blade which we interpreted as 
an initial position measurement with width 1 /ko/3. Damp- 
ing stems from the impurity-bath coupling which induces 
the kernel IT4l: 



TkiO) 



K\k\\Uk\^ cosh[u\k\{ph/2 - id)] 



2T:h 



fimh[u\k\pn/2] 



(21) 



with 9 = s — ir. In order to understand the effects induced 
by the non linear impurity-LL coupling in ([5]) we expand 
eq. (PP]) to second order in q. The result can be found in 
[12] and the correlation function can be calculated: 

Cit,t) ^ n'{t) - n*f3 C^<^{tr + ^ , (22) 

Here, Tl{t) and C°'^{t) arc the response and equilib- 
rium correlation functions in the high temperature limit 
hprij <ti 1 which prevails in the experiment. In the 
Laplace domain they read TZ{z) = (1/M£)[z^ -I- zu>cCe{z) + 
(17J)2]-i and C'"i{z) = (l//3z)[l/«;* - 7^(^)]. Linear re- 
sponse and correlator depend only on the "damping ker- 
nel" 



a{t) 



duj 



pMi 
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(23) 



As we pointed out, the final equilibrium value should be 
rather I/k/3 [12] than 1/k*/3 that would follow from ([22]). 
We conclude that, while the Gaussian approximation of 
([20]) [see [T^ for details] yields a realistic description of 
the impurity dynamics for short times, it cannot deliver 
the right correlation function for large times. However, it 
is straightfoward to construct an approximate correlator 



4 n* p 



+ (1- 



-Tnit\ 



1 1 

K/3 K*I3 



(24) 



which interpolates between the two asymptotic expres- 
sions ([^^ and p^ . Here, F is the effective damping 
induced by the Luttingcr bath. For small to moderate 
damping it is given by F ~ |/i(ri//a;c)e~^^^/'^'= [12]. We 
expect (|24p to be a realistic approximation of the impurity 
position width. 

Discussion. In [10] ^^K atoms play the role of the impu- 
rities moving in optical ID tubes through a Luttingcr liq- 
uid made of ^^Rb atoms. Both the ^^K and the ^""Rb are 
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confined in the same longitudinal optical potential with 
the (bare) potential frequency ilj ~ 550s^^ (390s~^) for 
^^K (^^Rb). We interpret the initial localization (with 
Kq ~ 150k) of the impurities as a position measurement. 
The experimental temperature is such that hftjf] ~ 10~^ 
which ensures the high temperature regime. The mean 
squared velocity is obtained to be « 8.5 mm/s which 
exceeds the typical sound velocity w w 3 mm/s so that the 
impurity moves in the supersonic regime. 

As pointed out before, the mass has to be renormalized 
in such a way that remains approximately constant 
over a wide range of w/wl- This can be achieved by a 
suitable choice of uJc, the only free parameter in our the- 
ory. For 7 = 0.35 the choice ojc/^i = 40 — 50 leads to a 
variation of 10 % for in the range < w/wl < 5. How- 
ever, since the mass and potential shifts decrease during 
equilibration the oscillation frequency can slightly change 
in time. Thus, for large times fij approaches flj in any 
case. Hence, in order to experimentally observe the n*j 
predicted by Eq. (|T6| one cannot average over many peri- 
ods as was done in [10] . It is therefore not straightforward 
to make a direct precise quantitative comparaison between 
Eq. (fTB|) and the experimental findings, although we think 
that the evidence in [T0| clearly indicates that the mass 
renormalization counteracts the potential shift to a large 
extent. 

Finally, we compare ([24]) to experimental data in Fig. [2] 
In [TU] \JC{t, t) has an offset of about 5/xm which wc add 
to our theoretical results. Moreover, for small interactions 
(w/wl ^1) there is a residual damping in the experiment 
due to inter-impurity collisions in tubes with several im- 
purity atoms [TO] which is of course not covered by our 
theory. We therefore use the data from [T^ for the damp- 
ing constant (F « 0.03 for w/wl ~ 1) in a{t). As pointed 
out before fl} can slightly vary during the equilibration 
process. However, this effect is not expected to be observ- 
able within the experimental error bars and therefore we 
approximate ^l} by ftj for all times. The match between 
the experimental data and our theoretic curves is quite 
impressive. 

In a recent paper jl2j we argued that non-trivial effects 
on ri/, produced by the super-Ohmic spectral density in 
the damping kernel (^5)) . could be observed for uJc ~ 0.3fi/. 
Now, after having gained a deeper insight into the impu- 
rity dynamics we think that Wc should be much larger such 
that the effects of the polaronic mass shift and the po- 
tential renormalization (which were previously neglected 
in [12] ) largely dominate the influence of the non-Ohmic 
spectral density. 

In summary we gained a thorough theoretical under- 
standing of the experimental data in [10] . We calculated 
the effective potential spring constant with an EOM ap- 
proach, which is expected to be correct for short times, 
and we obtained a result without any undetermined pa- 
rameter [see Fig. |T] . We argued that due to memory ef- 
fects neither a potential nor a mass renormalization can 
take place in the long time limit. One question not re- 



solved yet concerns the precise mechanism that links the 
mass and potential shifts which we hope will be revealed 
by future experiments. Finally, using the analytic results 
in }12) for the Brownian motion of a particle coupled to an 
exotic environment, after an initial position measurement, 
and with a phenomenological correction to the asymptotic 
limit, we described the experimental data for C(t,t) very 
accurately. 

* * * 
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